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Abstract 

The response of the small scales of isotropic turbulence to periodic large scale forcing is studied 
using two-point closures. The frequency response of the turbulent kinetic energy and dissipation 
rate, and the phase shifts between production, energy and dissipation are determined as functions 
of Reynolds number. It is observed that the amplitude and phase of the dissipation exhibit non- 
trivial frequency and Reynolds number dependence that reveals a filtering effect of the energy 
cascade. Perturbation analysis is applied to understand this behavior which is shown to depend on 
distant interactions between widely separated scales of motion. Finally, the extent to which finite 
dimensional models (standard two-equation models and various generalizations) can reproduce the 
observed behavior is discussed. 
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I. INTRODUCTION 


Statistical transients in turbulence remain a major challenge to both theory and modeling. 
The mechanisms by which a turbulent flow readjusts to new conditions, for example in 
boundary layers with sudden changes in wall roughness or pressure gradient 1 , are not entirely 
understood and continue to resist prediction by models. 

Another class of statistically time-dependent turbulent flows is defined by the presence of 
periodically oscillating forcing. The classic example is steady pipe flow with small superposed 
oscillations of the mean pressure gradient. This flow has been the subject of extensive 
experimental 2,3 theoretical 4,5 and numerical 6 investigation. There are two obvious limits: 
the ‘static’ limit of slow oscillations, in which the turbulence evolves through a sequence 
of local steady states, and a limit of ‘frozen’ turbulence in which the turbulence does not 
respond at all to the oscillations. 

Analysis of oscillating pipe flow typically concentrates on the phase relations among the 
wall shear, centerline velocity, and pressure perturbation. These quantities prove remarkably 
difficult to predict at frequencies intermediate between the static and frozen limits even if the 
problem admits a linearized description, indicating unanticipated subtleties in the dynamics; 
indeed, the only entirely adequate predictions are by Large Eddy Simulation 6 , which is very 
surprising in view of the apparent simplicity of the problem. 

Recently, the problem of periodically forced homogeneous isotropic turbulence has been 
proposed' and investigated theoretically , 8 by numerical simulations , 9 and by experiments 
using time-dependent grids 10 . Because of the absence of complications like near- wall be- 
havior, this problem provides an ideal setting in which to investigate the time-dependent 
spectral dynamics of turbulence. 

Previous work on this problem has been motivated by a search for resonance- like energy 
response near a critical frequency proportional to the inverse large-eddy turnover time, and 
perhaps at integer multiples of this frequency as well. This paper focuses instead on the 
properties of the dissipation rate. At frequencies intermediate between the static and frozen 
turbulence limits, nontrivial Reynolds number dependent properties are found. The energy 
cascade acts as a filter that suppresses oscillations at small scales, but as in the oscillating 
pipe flow, the details are more complex than the simple problem statement would suggest. 

The main results are obtained by the Eddy Damped Quasi-Normal Markovian (EDQNM) 
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closure 11 ' 12 . The predictions of this closure for periodically forced turbulence are in reason- 
able qualitative agreement with existing results. Elementary arguments show that at forcing 
frequency u, the amplitude of the energy and dissipation rate oscillations vary as cu _1 for 
large frequencies. But the calculations show that the dissipation rate modulation amplitude 
exhibits nontrivial scaling in the intermediate frequency range, and the phase difference 
between the production and the dissipation rate has complex dependence on both uj and 
Reynolds number in this range. 

To understand this behavior, we apply asymptotic analysis to two simpler models: the 
classical Heisenberg model 13,14 and a recent generalization 15 . In these models, the details of 
triad interactions are suppressed, but the essential idea of nonlocal interaction is retained. 
We show analytically how the energy cascade filters the oscillations, and that this filtering 
is responsible for the observations. 

Some finite dimensional models of the two-equation type will be considered. The two- 
equation model is correct in both the static and frozen limits, but misses important features 
of the dynamics at intermediate frequencies, including the Reynolds number dependence 
of the dissipation. A more complex three-equation model allows for more complex phase 
relations, but is also incapable of capturing the Reynolds number dependence. It should be 
emphasized that this dependence is not a low Reynolds number effect. 

A related problem involving periodic forcing is oscillatory homogeneous shear flow 16 . 
This problem has important, and even dominant linear effects; it therefore has a somewhat 
different character from periodically forced isotropic turbulence, in which only nonlinear 
mechanisms arc important. Another related problem can be mentioned, in which turbulence 
is forced periodically at the boundary of the flow region 17,18 . This flow has many interesting 
similarities to periodically forced isotropic turbulence; although it is simpler in many respects 
than periodically modulated pipe flow, the dynamics of this problem may include effects of 
turbulent diffusion as well as energy transfer and may therefore not be entirely amenable to 
the present type of analysis. 
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II. DEFINITIONS AND ELEMENTARY PROPERTIES 


The spectral evolution equation for time-dependent forced homogeneous isotropic turbu- 
lence is 14 

E(k, t ) = P(k, t ) — T(k, t ) — 2 uk 2 E(k, t ) (1) 

where E(k, t ) is the energy spectrum and T(k, t ) is the energy transfer due to nonlinear inter- 
actions. The production spectrum P(k, t ) is assumed to be localized near some wavenumber 
Kp(t). Consider a basic steady state, defined by the time-independent form of Eq. (1) 

0 = P(k) — T(k) — 2 uk 2 E{k). (2) 

The problem of periodically forced turbulence is formulated by introducing a periodic per- 
turbation of the production spectrum, 

P(k, t ) = P{k) + P(k) cos(u)t) (3) 

where we will assume 

P(k) = £P{k) (4) 

with e < 1, so that the problem can be analyzed by linearization about the steady state 
defined by Eq. (2). Then 

E(k, t) — E(k) + 8E(k, t) (5) 

with SE(k, t) <C E(k). If at sufficiently long times, E(k, t ) becomes periodic in time, linearity 
implies that the period is u, hence 

5E(k, t) = E(k) cos(u)t + 4>e{k)). (6) 

In terms of the quantities 

F(k) = E{k) cos (J)e(k-) G(k) = E(k) sin 4>e{k) (7) 

6E(n,t) is written as 

SE(K,t) = cos(c ot)F(n) — sin (ut)G(K). (8) 

The basic time dependent single-point moments: total production Pit), turbulent kinetic 
energy k(t), and dissipation rate e(t), are expressed in terms of their time averages P, k, 


4 



and e and their phase averages P, k, and e as 


where 


m 

= P + P cos (cat) 

( 9 ) 

k(t) 

= k + k cos(c ot + fa) 

(10) 

e(t) 

= e + e cos(c ot + fa) 

(11) 


P — dn P(k) P — / dn P(n) 

Jo Jo 

poo poo 

k — dn E(n) e — / d/c 2 un 2 E{n) 

Jo Jo 

and in view of Eq. (8), 

poo poo 

kcos(cj)k) = / dn F{n) ksin(fa) = / dn G(k) 


( 12 ) 


(13) 

(14) 


ecos(0 e ) = / dn2vn 2 F{n) esin(0 e ) = / dn 2 uk 2 G(k). 

Jo Jo 

For simplicity of notation, the spectral densities P(k) and P(k) are distinguished from the 
corresponding single-point moments P and P by their arguments rather than by a new 
letter. 

The simplest formulation of the problem seeks the dependence of the phase averaged 
amplitudes k and e and the phase shifts (f) k , 4>e on the forcing frequence cu; k will be called 
the modulated energy and e the modulated dissipation. Pit), k(t), and e(t) are related, 
independently of any closure hypothesis, by the energy balance, obtained by integrating Eq. 
(1) over all wavenumbers: 


k(t) = P(t) — e(t) 


(15) 


where energy conservation by nonlinear interactions implies that 

dn T(k, t ) = 0. 

Substituting Eqs. (9)— (11) in Eq. (15) and subtracting the steady balance P 
relation for modulated quantities 

— c ok sin(caf + fa) = P cos(c ot) — e cos(a >t + fa) 



(16) 
e gives the 

( 17 ) 
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or equivalently 


— uk sin <f>k = P — e cos <p e 
— c vk cos 4>k = e sin 

Elementary trigonometric identities give the explicit relations 


(18) 

(19) 


1 / 2 / ~ \ 2 
k — -W(e sin 4> e ) + (^P — e cos <p e 


tan (pk = 


P — e cos < 


e sm < 


( 20 ) 


and the equivalent relations 


e = \l [uk cos 

uk cos (pk 


2 / _ _ \ 2 
+ ( P + uks'mcpk 


tan® f = 


( 21 ) 


P + wfcsindfc 

Although additional assumptions are obviously required to close the problem, explicit 
closure hypotheses are not required to reach some simple but useful conclusions about the 
limits of asymptotically high and low oscillation frequency. Linearity implies that the fre- 
quency of the perturbation at any scale of motion must be the imposed frequency u, but 
in the inertial range, disturbances are damped on the Kolmogorov time-scale (e 1//3 K 2 / 3 ) -1 ; 
accordingly, we anticipate that if u 3> e 1//3 ft 2//3 , the perturbations must be over-damped, but 
that they are active and only weakly damped if u <C e 1//3 ft 2 / 3 . This argument suggests that 
in the static limit u | 0, the turbulence follows the slow modulations at all scales of motion, 
so that also (p e (u), (pk(u) | 0. Then Eq. (18) gives P ~ e; Eq. (19) is not satisfied exactly, 
but is approximately true since u ~ 0. Assuming that for slow modulations, the relation 
e(t) = C e k(t) 3 / 2 / L remains valid with time-independent L, and that the small perturbations 
k and e are nearly static, then e/e = (3/2 )k/k . These observations suggest that in this 
limit, the single-point modulated quantities admit series expansions in positive powers of u: 



(22) 

P + 0(u 2 ) 

(23) 

0(u) 

(24) 

0(u) 

(25) 
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where the powers of lo are suggested by the parity properties of Eqs. (18) and (19) under a 
change of sign of lo. Equivalently, to lowest order, we have 


e 


uk; 



(26) 


where the frequency Co defined by this equation is the ‘critical’ frequency discussed by Lohse 7 . 

In the ‘frozen turbulence’ limit lo \ oo, we see that Eq. (18) is satisfied if cpk ~ — 7 t /2; 
then k ~ P/io. If, as the simple argument above suggests, the perturbations are overclamped 
throughout the inertial range, the only scales of motion at which the oscillating force can 
be effective are the forcing scales themselves. If so, the modulated dissipation will also take 
place in this range of scales, so that 


e = 


dhi 2 vk 2 E(k) ~ 2 vk 2 p / dn E(k) ~ 2vn 2 P k ~ 2un, 


LO 


(27) 


As in the previous limit, Eq. (19) is not satisfied exactly, suggesting that the perturbation 
quantities should admit series expansions in negative powers of lo: 


k = - + O(io- A ) 

LO 


€ = 2 VK 


P 


P b 0{lO 

LO 


- 3 \ 


7 T 


4>k ~ + O(i0 X ) 


71 


= ~ +0^) 


(28) 

(29) 

(30) 

(31) 


III. NUMERICAL RESULTS FROM SPECTRAL CLOSURE 


In this section, we apply the EDQNM spectral closure 11 to this problem. The exact 
formulation of the model and the numerical method is the same as in Touil et al. 19 and 
for details we refer to that work. In this closure, nonlinear interactions among wavenum- 
ber triads of different ‘shapes’ are considered explicitly, with a definite weighting derived 
perturbatively from the governing equations. 

The energy spectrum was initialized by a von Karman spectrum; however, the influence of 
the initial energy spectrum vanishes after a transient and the results reported are evaluated 
after reaching an asymptotic state. The large scale forcing is 

P(n,t) = a (p + P cos(iot)^j exp(— 7 « 2 ), k > k 0 > 0 (32) 
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FIG. 1: top: A; as a function of lo/uj for R\ varying from 30 to 1000; bottom: same for cok 

/»oo 

where a is the normalization constant such that a / exp(— 7 K 2 ) = 1, with P/P = 0.125 

J K 0 

and 7 = 0.5. The spectral resolution is approximately 20 wavenumbers per decade. The 
results are shown in Figures 1-4. 

1. Modulated kinetic energy k. 

Figure 1 shows a plateau for k at low frequencies and cu -1 dependence for high frequencies 
as suggested by the elementary arguments leading to Eqs. (22) and (28): the static and 
frozen turbulence limits are well reproduced. In our calculations, no local maximum of k is 
observed. In the DNS results of Kuczaj et al. 9 a small local maximum was present around 
the turbulent frequency u defined in Eq. (26), but in the shell model study by von der Heydt, 
Grossman and Lohse 20 , this maximum was absent. The existence and explanation of this 
maximum remain open questions. However, all of the available data exhibits a clear response 
maximum of the compensated quantity c ok near G). This maximum is also prominent in the 
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FIG. 2: top: e as a function of uj/uj for R\ varying from 30 to 1000; bottom: same for uie. 

EDQNM results shown in Figure 1. We leave the question of whether a response maximum 
of k itself is or is not consistent with closure unanswered for now. Conceivably, the answer 
is not universal, but may depend on the forcing scheme. The Reynolds number, or viscosity, 
does not seem to play an important role for k\ for moderate and high Reynolds numbers, 
all the data collapses on a single curve. 


2. Modulated dissipation e. 

Figure 2 shows that e also displays a plateau in the static limit as predicted by Eq. (23). 
Like the compensated quantity uk, the compensated data cue shows a response maximum 
approximately near uj. Beyond this frequency, e decreases sharply; at high Reynolds number, 
e ~ u~ 3 . But at even higher frequencies the cn _1 frequency dependence predicted in Eq. 
(29) is observed. It is interesting to note that the high frequency u; -1 range depends on the 
Reynolds number, and is indeed proportional to the viscosity, as suggested in Eq. (29). 
What remains to be explained is the fast drop of e at intermediate frequencies. Intu- 
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FIG. 3: top: Phase lags —(f>k as a function of uj/u for R\ varying from 30 to 1000; bottom: Phase 
lags — 4> t for different R\ 



FIG. 4: — (cj) e — 4>k ) as a function of lu/lu for R\ varying from 30 to 1000 

itively it can be explained as follows: at low frequencies the energy cascade can follow the 
modulation. At high frequencies the cascade filters the modulated energy flux, since the the 
turbulent frequency is lower than the modulated frequency. The fast drop corresponds to 
the rate at which the energy cascade filters the modulated energy flux. Insights into this 
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process have important physical consequences as they clarify how small scales are influenced 
by large scale forcing. 

3. Phase shifts. 

Phase shift data is shown in Figures 3 and 4. The phase lags <f>k and <f> e both go to zero 
for small uj. In this limit, everything is in phase as suggested by Eqs. (24) and (25). At 
high uj, (fk and <f> e go to 90 degrees, consistently with Eqs. (30) and (31). A slight overshoot 
in 0^ is observed around the turnover frequency. At intermediate values <f> e shows a large 
overshoot with respect to 90 degrees and a very noticeable dependence on Reynolds number. 
This can be explained as follows: as long as the energy cascade can follow the modulation, 
i.e. at low frequencies, the modulated energy is transfered to the dissipation range through 
the energy cascade. The finite cascade time T c introduces a phase shift between k and e 
proportional to ujT c . This is illustrated in Figure 4. At low frequencies <p e — (pk is a linear 
function of uj, which permits determining the cascade time. 


IV. ANALYTICAL TREATMENT BY SPECTRAL CLOSURE 

We supplement these numerical computations with analytical results. The complexity 
of the EDQNM transfer integral does not permit simple direct analysis, so we will consider 
much simpler models which embody certain features of nonlinear turbulence dynamics, but 
in a way that permits analytical conclusions to be drawn relatively easily. 

A. General formulation 

The general closure equation is found by introducing the closure hypothesis 

T( K ,t) = (33) 

in Eq. (1). Eq. (33) expresses the energy transfer in terms of the energy flux fF , which is 
assumed to be a functional of the energy spectrum. In the problem of periodic forcing, the 
perturbation 6E(n,t) defined by Eq. (5) satisfies 

5E(k , t) = P(k) cos(ut) — C[5E(n, t)] — 2 vk 2 5E(k, t ) (34) 
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where C is the linear functional 

m.,t)} = ^(^pj_m K , t )] ( 35 ) 

and (5T /8E)e denotes linearization of T at the steady state E(k). 

Separating terms proportional to cos(u >t) and sin(u;f), Eq. (34) can be written as 


■uG{k) = P{k) - C[F(hi)} - 2 vhi 2 F(hi) 

(36) 

■uF(k) = jC[G(hi)] + 2uhi 2 G(hi). 

(37) 


In view of Eq. (35), 

dhi £[<h(ft)] = 0 (38) 

therefore integration of Eqs. (37) and (36) recovers the single-point relations Eqs. (18) and 
(19). 

Before beginning the analysis, we note that substituting Eqs. (36) and (37) in Eq. (13) 
gives 



1 _ i 

k sin(0fc) = PH — / 2vK 2 F(n)dK 

u u J 0 

1 r°° 

kcos(cj)k) = / 2 un 2 G(n)dn. (39) 

Jo 


Ignoring the viscous terms recovers k sin 0 fc « —u~ 1 P, which is equivalent to k ~ lu~ 1 P and 
( f> k ~ — 7t/ 2, the approximations obtained by elementary arguments as Eqs. (28) and (30). 
The corresponding substitutions in Eq. (14) yield 


ecos®, = — - 


e sm 0. = — - 


U) 

1 

OJ 


Uo 


00 


dhi 2 vk 2 C[G{k)\ + dhi 4 v 2 hi 4 £.[G(hi) 


dhi 2 vhi~P(hi) — dhi 2 uk, 2 C[F(hi)\ — dhi 4v 2 hi 4 jC[F ( n)} 


(40) 


.(41) 


Obviously, very strong assumptions are needed to reach any conclusion about e and 0 £ , 
demonstrating that the behavior of the oscillating dissipation rate is somewhat subtle. Thus, 
the elementary conclusion that e can be approximated by taking only the first term in Eq. 
(41) requires arguing that the terms in v£, which represent oscillatory vortex stretching, 
can be ignored, and that, despite the presence of n 4 in the corresponding integrals, the 
terms in z/ 2 , which represent oscillatory enstrophy destruction, can also be neglected. These 
assumptions are much less convincing than those underlying the elementary approximation 


12 



for k. In fact, more careful analysis will reveal nontrivial features of the dynamics of the 
modulated dissipation. But these features can only be computed using a model; this issue 
will be considered in the next section. 

B. Simplified integral closure models 

Kraichnan 21 showed that if the correlation equation in a closure of the DIA family is 
simplified by restricting attention to distant interactions only, an energy transfer model 
close in structure to the classical Heisenberg model is obtained. We refer to 21 for detailed 
explanations. Following this observation, Rubinstein and Clark 15 constructed a generalized 
Heisenberg model by adding asymptotically local interactions to the transfer model. The 
result is the energy flux closure 

E[E(k)\ — C | J dp p 2 E(p) J dp E(p)9(p) — J dp p 4 J dp — — - — j> . (42) 

where the time argument is not explicitly written. Here and subsequently C will denote 
some constant, but not necessarily the same constant each time it appears. This energy 
transfer model was supplemented 15 by an evolution equation for the time-scale 9(k), but the 
present work will use the simple algebraic closure, 

9(E) = [k 3 E(k)]~ 1/2 . (43) 

Theoretical features of this model include the possibility of energy transfer from small to 
large scales (energy ‘backscatter’) due to the possibility of negative energy flux, and con- 
sistency with the existence of inviscid equipartition ensembles in which E(k) oc k 2 . We 
note that Canuto and Dubovikov 22 had already obtained a simple spectral model consistent 
with many of the properties of analytically much more complex models by supplementing a 
renormalization group model with a backscatter term. 

The classical Heisenberg model 13 is obtained by discarding the negative term in Eq. (42), 
so that 

PK, POO 

E[E{k)\ — C / dp p 2 E(p) / dp E(p)9(p). (44) 

J 0 Jk, 

In this model, the energy flux is necessarily positive; hence energy is always transferred 
from large scales to small scales. This property is inconsistent with the possibility of in- 
viscid equipartition. Despite these drawbacks, the Heisenberg closure models one feature 
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of turbulent energy transfer that will be crucial to the present analysis: the possibility of 
‘distant’ interactions between modes with disparate wavenumbers. 

The linearized transfer for the generalized Heisenberg model is 


£[$(«)] = C < k 2 $(k) / dp 


E{p) ' [E{d 


-k 2 E(k) / dp 


$(p) 


1 5>(/c) 


{ E(p)p 3 y / 2 2{E(k)k 3 } 1 / 2 J 0 


dp p 2 $(p)+ 


dp p 2 E(p) 


-k 4 / dp 


E(p) l / 2 $(p) 3h(fi) 1 / 2 $(fv) 


P 


,7/2 


+ 


/H/2 


dp p 4 


and the linearized transfer for the classical Heisenberg model is 

poo 

£[$(«)] = C ^ /7 2< h(/7) / dp 

J K, 

$(p) 


E(p) / E(k) 


-k 2 E(k) / dp 


po y K° 

1 <!>(«;) 


dp p 2 $(p)+ 


dp p 2 E(p) 


(45) 


(46) 


{E(p)p 3 } 1 / 2 2 {E^K?} 1 / 2 

Returning to the analysis of Eqs. (36) and (37), we note that they can be decoupled to 
give 


(u 2 1+(jC + 2uk 2 1) 2 )F(k) = jC[P(k)} 

(cu 2 l+ (C + 2 uk 2 1) 2 ) G(k) = -ujP{k) (47) 

so that inversion of the linear operators on the left gives the solution for F and G. But 
since exact inversion is only possible numerically, we will seek asymptotic solutions for large 
u using standard methods. A lowest order approximate solution of Eqs. (36) and (37) is 
obtained by balancing the leading order terms in ui, so that 

F(k) « 0 G(k) « -P(k). (48) 

U) 

Since this result ignores nonlinearity, it might be called ‘rapid distortion theory’ for this 
problem. 

A formal solution of Eqs. (36) and (37) can be constructed in powers of by perturbing 
about the leading order solution Eq. (48); taking only the correction terms of the next order 
gives 


F(k) = c o~ 2 C[P(k)} 

G(k) = -u~ 1 P(k) +uj- z P 2 [P(k)\. (49) 
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This approximation can also be obtained by operator inversion in Eq. (47) by a Neumann 
series. The resulting series is divergent but asymptotic in therefore, as usual in such 
cases, the truncated series can provide useful information. 

The corrections in Eq. (49) depend on C[P(k)]. We note from Eqs. (46) and (45) that all 
contributions to jC[P(k)} are proportional to P(k), and therefore vanish when -P(ft) vanishes, 

/ OO 

P{jp)dp and vanish for large n since P(k) is nonzero only for small 


J K, 

k, with one exception, the term common to both models, 


C nl [P(k)] = i/S f dp p 2 P(p), (50) 

v K Jo 

where the subscript NL denotes that these are contributions from nonlocal interactions. 
Because this term pertains to forward transfer alone, we obtain it in the Heisenberg model; 
the backscatter term in the generalized Heisenberg model therefore plays no role in this 
particular analysis. 

The contribution to linearized transfer in Eq. (50) shows that the oscillatory disturbance 
is not confined to the region where P(k) is nonzero, even at asymptotically large u, contrary 
to the conclusion suggested by elementary considerations. Instead, the oscillatory distur- 
bance can propagate into all scales of motion. The remaining terms in Eqs. (45) and (46) 
with <3 >(k) replaced by P(k) can be considered corrections to the leading order solution Eq. 
(48), and will be ignored in the following analysis; the term in Eq. (50) itself provides the 
leading order solution in the regions where P(k) vanishes. 

For large k, 

C nl [P(k )] = J - [ dp p 2 P(p ) ~ KpPe 1/3 «T 7/3 
" K Jo 

C 2 nl [P(k)] ~ \f^£ dP K 2 P Pe 1/3 p- 7/3 ~ 4Pe 2/3 R- 5/3 . (51) 

Adding these nonlocal contributions to the leading order solution Eq. (48), the approxima- 
tion Eq. (49) takes the form 


F(k) ~ U- 2 K 2 p Pt l3 K~ 113 

G(k) ~ w- 1 P(K)+w- 3 4Pe 2/3 P 5/3 . (52) 

Eq. (51) shows that C p nl [P(k)] ~ k L0 ^ 3+p ] it follows that the series expansions for F(n) 
and G(k) proceed in positive powers of K] the higher order approximations will eventually 
contain positive powers of k. indicating the divergence of the series noted earlier. 
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Assuming that scaling ranges for F and G both begin at a scale of the order of n P , we 
will have 


k cos(0k) = / d,K F(k) ~ lu 2 e 1 ^ 3 n p 3 P 


k sm(4>k) — dn G(k) ~ — lu l P + lu 3 e 2 ^ 3 n 4 p 3 P. (53) 

Jo 

If lu is large, the terms of order larger than u;” 1 can be ignored, and we again return to the 
elementary estimate k ~ P /lu and (f)k ~ 7r/2 with corrections depending on /lu. 

The situation is quite different for the modulated dissipation rate, for which 


ecos (4>k) = 
esin(0 fc ) = 


dn 2uk 2 F(k) ~ lu 2 PJpPP^vKp = lu 2 P 2 p PP^ 2 v v ^ 2 


dn 2 vn 2 G(n) ~ lu 1 2uk p P + lu 3 n p Pe (54) 

•to 

where Kd = (e/z/ 3 ) 1 / 4 is the Kolmogorov scale. Evidently, there is a competition between 
the limits lu — > oo and v ~ Re~ 3 — > 0. The limit lu — > oo at fixed Re will indeed recover 
the elementary result e ~ o> -1 , but at hxed large lu, the limit Re —> oo gives instead e ~ 
LU~ 3 n 2 P Pe. The phase has the general approximate value tan (j) e ps c uv l ^e~ 1 ^ 2 + lu^ 1 u~ 1 ^ 2 P^ 2 
indicating a complex joint dependence on lu and Re in general. 

The main consequences of this analysis are the cu -3 range for e and the complex depen- 
dence of </> e on Reynolds numbers; both are confirmed by the EDQNM calculations. 


C. Scaling analysis for e 


The analysis in the previous section shows how nonlocal interactions in the Heisenberg 
and generalized Heisenberg models can carry the oscillatory disturbance into the inertial 
range. These observations suggest a simple scaling analysis for the modulated energy flux. 
Assume, following the discussion in Sect. 11 that scales of motion for which the oscillations 
are overdamped, that is, scales satisfying 9{n) 1 > lu do not transfer any modulated flux, but 
that modulated flux is transferred by scales of motion such that 9{np l < lu. The crossover 
occurs at the scale n u defined by 9(n)~ l = lu, or n^ = In both the Heisenberg and 

generalized Heisenberg models, the transfer of modulated flux is then given approximately 
by 


r^p 


dp p 2 E(p) / dp E{p)d{p) ~ ttpfce 1 / 3 *^ 4 / 3 


kIPlu-^lu - 2 ?' 3 


k%,Plu 


- 3 e. 


(55) 
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This result is consistent with the existence of a contribution to G scaling as k - 5 / 3 obtained 
more formally in Eq. (52). 

The argument can be extended to the EDQNM closure as follows. Modulated kinetic 
energy is injected in the flow around the wavenumber Kp. This energy will leave the large 
scales to enter the energy cascade at a rate e m (ftp). Using classical reasoning, this rate can 
be estimated by: 


e in {n P ) ~ 


k 

9(k p ) 


(56) 


at high frequencies the modulated energy is: 


k ~ Pu~ l 


(57) 


and the timescale can be estimated by 

9 (up) ~ e _ 1 / 3 Kp 2 ^ 3 


(58) 


so that 

e in (K P ) ~ Pa; _1 e 1/3 Kp 3 . (59) 

The point is now that this energy will be overdamped if it passes through the scales Kp < 
k < k u . The only way to reach the zone that can transfer the modulated flux, k > k^ is by 
nonlocal energy transfer. This transfer will involve, for Kp <C triads with two legs of a 
length Ku and one leg equal to Kp. The disparity parameter s defined as 


max(K,p, q) 
min(/c,p, q) 

with k, p, q the norms of the wavevectors forming a triad, is for these triads 

Kuj U 3/2 

Q ^ ~ 

Kp Kp6 1 / 2 


(60) 


(61) 


It was predicted by Kraichnan 23 (compare also the DNS study by Zhou 24 ), that the nonlocal 
part of the energy transfer involving triads with a disparity around s, €-'(k,s) with respect 
to the total energy flux €^(k) scales as: 



; - 4/3 


(62) 


In our case we identify the total flux of modulated energy, epK,s) with e m (Kp). The 
nonlocal flux e*(K,s) corresponds to the modulated energy flux that manages to reach the 
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range k > and that will eventually be dissipated, and so is equal to e. One finds therefore 
combining Eqs. (56) and (62) that the modulated dissipation for high frequencies equals: 

e ~ ~~~ W 4 / 3 ~ K 2 P Peuj~ 3 (63) 

0{Kp) 

in agreement with Eq. (55). The inviscid nature of this correction to the modulated dissi- 
pation is in agreement with the observation in Figure 2. 

An important distinction between the classical and generalized Heisenberg models and 
EDQNM is that the power-law scaling of Eq. (62) applies for all s in the simple models, but 
is given by a more complex expression for EDQNM. This implies a difference in the detailed 
predictions when is of order one. 


V. FINITE DIMENSIONAL MODELS 

The problem of periodically forced turbulence has been investigated through properties 
of the single point moments k(t) and e(f); complete results for these quantities have been 
found from various spectral closure theories. Single-point modeling attempts to circumvent 
spectral modeling by constructing closed equations for the single point moments themselves. 
It is an important theoretical question whether such equations exist 25 , and indeed, much 
stronger assumptions are needed to close the problem at this level. In this section, we will 
assess how much of the dynamics is accessible to single-point modeling. 

In order to permit the underlying steady state, a two-equation model for periodically 
forced turbulence must take the form 


k = P — e (64) 

e = cUP-e) (65) 

k 

where Eq. (64) is just the energy equation previously stated as Eq. (15). For forcing at a 
fixed length scale, it can be shown 26 that C = 3/2; the e transport equation Eq. (65) then 
states that L = k 3 ^ 2 /e is constant, since Eqs. (64) and (65) imply L/L = (3/2 )k/k — e/e = 0, 
the same argument that gave Eq. (22). 

Eqs. (64) and (65) admit a steady solution in which P(t) — P — e(t) = e. We consider 
the perturbation about this steady state due to oscillating forcing Eq. (9); linearization 
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about the steady state and using the value of the model constant C — 3/2 gives 


— c ok sin (cut + 0 fc ) — P cos (cut) — e cos (cut + 0 e ) 
—cue sin(cut + 0 e ) = cu P cos(cut) — e cos(cut + 0 e ) 

where Eq. (66) restates Eq. (17). Divide Eqs. (66) and (67) to obtain 

sin (cut + 0 £ ) e 


so that 


and 


sin (cut + 4> k ) k 

0fc = 0e = ' 

e = uk 


= cu 


The linearized equations reduce to 


— cu/c sin 0 = P — uk cos 0 
— cu cos 0 = cu sin q b. 


( 66 ) 

(67) 

( 68 ) 

(69) 

(70) 

(71) 


Note that this is just the general result of Eqs. (18)-(19) with the special closure hypothesis 
4>k = 0 e . Indeed, this is perhaps even a rather plausible closure, since 0*, = 0 e is certainly 
true in both limits cu { 0 and cu | oo. It follows that 


tan (f> 



cu 


P P 1 

— cos <p = — — =. 
u u y/l + (cu/cu) 2 


(72) 


The limits 


0 ~ — 7t/2; k ~ P/cu for cu — » oo 

. _ (73) 

0 ~ 0; k ~ P/cu for cu — >• 0 

are consistent with the limiting results previously obtained as Eqs. (28) and (22). Whereas it 
is certainly expected that a two-equation model should be adequate in the static limit, it may 
be surprising that the frozen turbulence limit for k is also predicted correctly, particularly 
in view of the suggestion 5 that in oscillating channel flow, predicting the frozen turbulence 
limit requires rapid distortion theory. 

Despite these successes, the two-equation model has some important limitations. First 
of all, the phase shifts 0*, and 0 e are equal for all cu, in disagreement with Figure 3. Second, 
Eq. (70) states that e and k are proportional for all cu. Recall that this proportionality was 
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found in Eq. (26) as a consequence of assuming a constant length scale. Comparison of 
Figures 1 and 2 shows that k and e are certainly not always proportional. This comparison 
demonstrates that the identification of the constant forcing scale Kp 1 with a multiple of the 
ratio k 3 /’ 2 /e cannot be made for general values of u ;; following the common terminology that 
‘equilibrium’ turbulence is turbulence in which all dimensional arguments are valid, we can 
say that periodically forced turbulence is not in ‘equilibrium.’ 

What is most striking is that the two-equation model cannot predict the //-dependence of e 
and (f> € , which is not a low Reynolds number effect in this case. A fundamental observation of 
Speziale and Bernard 2 ' is that Reynolds number dependence in the dissipation rate dynamics 
is a manifestation of unbalanced vortex stretching , the absence of which underlies the classic 
formulation of the e equation by Tennekes and Lumley 28 . Even if we were satisfied with a 
high Reynolds number model, it should predict e ~ 0 for large u. We have noted that this 
limit is due to the filtering effect of the spectral cascade. Evidently, this effect cannot be 
captured at the level of a two-equation model. 

Another way to understand the relation between spectral closure and the two-equation 
model is to note that Eq. (71) is obtained from the general closure model Eqs. (36) and 
(37) by making the single relaxation time approximation 

C + 2uk 2 1 « ul. (74) 

before integrating over k. This type of simplification, by which the continuum of time scales 
in a turbulent flow is replaced by a single dominant time scale, is a mainstay of modeling, and 
is often very useful; however, in the problem of periodically forced turbulence, it suppresses 
the nontrivial features of the finite u dynamics. 

One remedy is, as always, to argue that the model constants should be functions. If we 
set C = C{u/uj), then if C | 0 for uj | oo, the correct behavior can be reproduced. However, 
this ad hoc model would have no validity apart from this very special problem and would 
merely amount to curve-fitting. 

We would like to comment briefly on the modeling of this flow with a more complex finite 
dimensional model with two characteristic time scales; that is, a ‘multiple-scale’ model 26 . 
For example, consider a three-equation model in which energy flux / is distinguished from 
dissipation e. A general form for such a model that is consistent with a steady state is 
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k = P-e 
f = c J-(P - f ) 

i = C 2 |(/- 1). (75) 

The limits C'i f oo and C 2 ] oo both recover the two-equation model. The primary mo- 
tivation for this model is that there are now two time scales /c/e and k/f instead of only 
one. 

Linearizing as usual about the steady state 


k{t ) = k + k cos(c ot + 4>k) 
fit ) = f + f COs((Vt + (pf) 

e{t) = e + e cos(u;f + 0 e ). (76) 

Then 


—uj k sin [ut + 4>k) = P cos (cat) — e cos (cat + <f > e ) 

I \~ 

— ca/sin(cat + 0/) = CW Pcos(cat) — /cos(cat + 0/) 

k L 

— caesin(cat + 0 e ) = CW /cos(cat + 0/) — ecos(c at + 0 e ) 

rC L 


(77) 

(78) 

(79) 


The intervention of the new quantity / in the dynamics means that k and e are no longer 
constrained to be in phase. However, prediction of the high Reynolds number result e ~ 0 
remains impossible: Eq. (79) then requires / = 0, which is inconsistent with Eq. (78). 
It is not difficult to evaluate both phase lags 0/ and 0 e , but even without explicit results, 
it is evident that Reynolds number dependence of 0 e remains inaccessible to this model. 
Although the three-equation model allows more complex phase relations and modeling of 
time delays in the spectral cascade, it, like the two-equation model, cannot take the Reynolds 
number dependence into account correctly. The addition of time scales to the dissipation 
rate dynamics does not solve all of the problems of two-equation modeling. 


VI. CONCLUSIONS 

The influence of periodic large scale forcing on isotropic turbulence was investigated by 
spectral closure theory. The asymptotic frequency dependence of the modulated energy and 
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modulated dissipation as observed in recent simulations 9 were recovered. It was pointed out 
that the asymptotic behavior of the modulated dissipation, which is proportional to cu _1 , cor- 
responds to the viscous damping of the forced wavenumbers, which is local in wavenumber 
space. For high and moderate Reynolds numbers an intermediate cu -3 frequency depen- 
dence of the modulated dissipation was observed in the EDQNM calculations. This range 
characterizes the filtering properties of the energy cascade. Closures allowing for nonlocal 
interactions (EDQNM, classical or generalized Heisenberg) can reproduce this behavior as it 
corresponds to nonlocal energy transfer between the forced scales and a range of wavenum- 
bers characterized by a crossover wavenumber k > ~ y / cu 3 /e. Finally it was argued that 

finite dimensional models can not correctly describe the problem of modulated turbulence. 
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